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Abstract. Fractal dimension is central to understanding dynamical processes occurring on networks; how- 
ever, the relation between fractal dimension and random walks on fractal scale- free networks has been rarely 
addressed, despite the fact that such networks are ubiquitous in real-life world. In this paper, we study 
the trapping problem on two families of networks. The first is deterministic, often called (a;, y)-flowers; the 
other is random, which is a combination of (l,3)-flower and (2,4)-flower and thus called hybrid networks. 
The two network families display rich behavior as observed in various real systems, as well as some unique 
topological properties not shared by other networks. We derive analytically the average trapping time for 
random walks on both the {x, y)-flowers and the hybrid networks with an immobile trap positioned at an 
initial node, i.e., a hub node with the highest degree in the networks. Based on these analytical formulae, 
we show how the average trapping time scales with the network size. Comparing the obtained results, we 
further uncover that fractal dimension plays a decisive role in the behavior of average trapping time on 
fractal scale-free networks, i.e., the average trapping time decreases with an increasing fractal dimension. 

PACS. 05.40.Fb Random walks and Levy flights - 89. 75. He Networks and genealogical trees - 05.60.Cd 
Classical transport - 89. 75. Da Systems obeying scaling laws 



1 introduction 

' Complex networks have become a powerful and common 
] tool for studying complex networked systems in nature 
and society, which allow for describing quantitatively the 
' structural features and complexity of real systems [TJ[2] . A 
, fundamental problem in the field of complex networks is 

■ to unveil the influence of topology on diverse dynamical 
' processes taking places on networks [5] ■ As an important 
. stochastic process, random walks are a subject of a large 
' volume of research [4l[5l|6l[7l|8] , due to their wide range 
I of applications in various areas [51 [TUllll| . A primarily in- 

■ teresting quantity related to random walks is the first- 
; passage time (FPT) pifnUIHITCl lTB ) . Recently, the mean 
. first-passage time (MFPT) to a target node on a network 

■ has received an increasing attention |17lll8lll9ll2"0ll21ll22l 
I [^[21] , since it can serve as a measure of search efficiency 

■ to find the target. 

A remarkable properties unveiled by extensive empiri- 
cal research is that many, perhaps most real-life networks 
display scale-free phenomenon, meaning that their degree 
distribution P{k) follows a power law P{k) ~ 
This property constitutes our basic understanding of the 
structural organization of real systems and has a profound 
impact on random walks occurring on scale-free networks. 
A lot of work showed that search efficiency is substantially 
improved if the target is a node with the highest degree in 
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scale- free networks [26ll27ll28ll29ll30ll31ll32ll33) . Neverthe- 
less, scale-free behavior cannot reflect all the structural 
information of real networks. It was acknowledged that in 
addition the power-law property, a variety of real-life sys- 
tems exhibit fractal scaling, which is often characterized 
by the fractal dimension [31]. Taking into account fractal 
scaling of scale-free networks can lead to a better under- 
standing of how the underlying systems work |i35j . It was 
shown that fractality is unfavorable for the target search 
problem in fractal scale-free networks [551137] . However, 
what is the relation between the fractal dimension and 
scaling of MFPT for random walks on scale-free networks 
remains not well understood. 

In this paper, we study the classic trapping problem 
on two families of networks, which is a random- walk is- 
sue with an immobile trap located at a fixed position ab- 
sorbing all walkers that visit it. The first network fam- 
ily is of a deterministic type, called {x, y)-flowers [55115^ 
that contains the intensively studied pseudofractal scale- 
free web pHiT] and fractal hierarchical lattices [I^I431I44[ 
I45j as their limiting cases by adjusting the parameters x 
and y. The second family is random that can be looked 
on as a mixture of (1, 3)-flower and (2,4)-flower and thus 
called hybrid networks ^IB] . Both network families present 
some typical topological properties as observed in real-life 
networks. We focus on a particular case of the trapping 
problem where the trap is fixed at an initial node (namely, 
node with the largest degree). We derive analytically the 
average trapping time for both the {x, j/)-flowers and the 
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hybrid networks, based on which we obtain the depen- 
dence relation of average trapping time on the network 
size. Our results show that in both network families, the 
average trapping time grows algebraically with their size. 
Lastly, we present a detailed analysis of the obtained re- 
sults, and show that the increasing fractal dimensions of 
fractal scale-free networks may lead to enhancement of 
search efficiency for those target nodes with the highest 
degree. 



2 Network models and their properties 

Here we introduce two families of networks and their struc- 
tural features. The first one is deterministic, while the sec- 
ond one is stochastic. 



2.1 Construction and topology of (a;, ?/)-f lowers 




Fig. 2. (Color online) Illustration for the evolution process of 
the (1, 3)-flower. 




Fig. 1. (Color online) Construction approach of the {x,y)- 
flowers. To produce the networks at iteration n + 1, we can 
replace each edge in the networks at iteration n by two parallel 
paths with lengths x {x > 1) and y (y > x and j/ > 1) on the 
right-hand side of the arrow. When x = 1, any pair of old 
nodes directly linked to each other by an old edge creates y—1 
new nodes. These y — 1 new nodes and the two old ones form 
a red path of length y\ while the old edge is remained. When 
X > 1, each old edge connecting two old nodes is deleted and 
simultaneously replaced by two paths with the two old nodes 
as the endpoints of both paths: the y—1 red nodes and the 
two old nodes form a path of length y, while the x — 1 green 
nodes and the two old nodes form the other path with length 



The first family of studied networks, frequently called 
(x, y)-flowers |38[I39) . are constructed in an iterative way. 
Let Fn{x, y) {n > 0) denote the (x, y)-flowers after n iter- 
ations. Without loss of generality we assume that x < y 
and y > 1. The (x, t/)-flowers are built in the following 
way, see Fig. [TJ For n — 0, Fo{x, y) consists of two initial 
nodes connected by an edge. For n > 1, Fn{x,y) is gener- 
ated from Fn-i{x,y) through replacing each exiting edge 
in Fn^i{x,y) by two parallel paths consisting of x and y 
links. Figures [2] and [3] illustrate the iterative processes of 
two special networks: (1, 3)-flower and (2, 2)-flower. 

The iterative construction allows for determining rel- 
evant properties of the (x, 2/)-flowers. It is easy to derive 
that the numbers of nodes and edges in Fn{x,y) are 



x + y 




Fig. 3. (Color online) Iterative growth for the (2, 2)-flower. 
and 

M„ = {x + yT , (2) 
respectively. And the average degree of Fn{x, y) is 



{k)r 



2M„ ^ 2{x + y - l){x + y) 
iV„ [x + y - 2){x + yy^-^ + I 



(3) 



y-1 



y-1 



which approaches 2{x + y — ^)/{x -f- y — 2) as n — )■ oo. 

The {x, y)-flowers display rich behavior in their topo- 
logical structure |38ll5U] . They follow a power- law degree 
distribution P{k) ~ k~'^ with the exponent 7 = H-ln(a:-|- 
?;)/ln2 belonging to the interval [1 -I- In 3/ In 2, 00). For 
X — 1, the (x, y)-flowers are small- world [37] but non- 
fractal; while for a; > 1, the networks are "large- world" 
and fractal with the fractal dimension df = hi{x + y)/ Inx. 
In addition, in the whole family of Fn{x,y), all networks 
with the same parameter x + y have the same identical 
degree sequence, hence the same degree distribution. 
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2.2 Generation and architecture of hybrid networks 



Different from the first family of networks, the second 
family of networks in question are built iteratively but 
randomly ^46j, see Fig. 21 Since the networks are a ran- 
dom mixture of (1, 3)-flower and (2, 2)-flower, we call them 
hybrid networks. Let Hn {n > 0) represent the hybrid 
networks after n iterations. Then the networks are con- 
structed as follows. For n = 0, the network Hq is composed 
of two initial nodes connected by an edge. For n > 1, Hn is 
derived from Hn-i- replace each link existing in Hn~i ei- 
ther by a square on the top right of Fig. |4] with probability 
g, or by a diamond on the bottom right with probability 
1 — q. For q = 1 and <? = 0, both networks are determinis- 
tic, which are separately (1, 3)-flower and (2, 2)-flower, see 
Figs. Hand El 




Fig. 4. Iterative construction method of the hybrid networks. 
Each edge is replaced by either of the connected clusters on 
the right-hand side of arrows with a certain probability, where 
black squares represent new nodes. 



3 Random walks with a trap fixed on an 
initial node 

Having introducing the networks, we investigate the sim- 
ple discrete-time random walks on Fn{x,y) and Hn- At 
each time step, the walker moves uniformly from its cur- 
rent location to one of its neighboring nodes. Here we con- 
centrate on a special random- walk issue with a trap placed 
at an initial node denoted by ix, which is also known as 
the trapping problem. One of the most significant quantity 
characterizing the trapping problem is the trapping time, 
i.e., the so-called first-passage time, which is defined as the 
expected time spent by a walker, starting from a source 
node, to arrive at the trap for the first time. The average 
trapping time, namely, mean first-passage time to the trap 
node, defined as the average of trapping time over all start- 
ing points other than the trap, measures the efficiency of 
the trapping process: the less the average trapping time, 
the higher the efficiency, and vice versa. 

Let rj*-"-* stand for the trapping time for a walker start- 
ing off from node i to first visit the trap ix in Fn{x, y) or 
Hn- Then, the average trapping time, (T)„, which is the 

mean of T-^"^ over all non-trap starting nodes in Fn{x,y) 
or Hn, is given by 

{T)n = ^ E , (4) 

where f2n denote the set of nodes in Fn{x, y) or Hn- It is 
obvious that for all n we have T*-"-' = 0. In the following, 
we will find the exact solution to (T}„ for both Fn{x,y) 
or Hn, as well as the dependence relation of (T)„ on the 
network size 7V„. 



The number of edges and nodes in Hn are M„ = 4" 
Nn = 2(4" + 2)/3, respectively. The average degree of 
Hn is {k)n = 2Mn/Nn = (3 X 4")/(4" + 2), which ap- 
proaches 3 in the limit of large n. In the full range of 
< g < 1, the whole family of networks has identical de- 
gree sequence and degree distribution with a power-law 
form P{k) ~ k~^- The clustering coefficient of the net- 
works is zero, since they contain no triangles. For q = 0, 
the network is "large-world" and fractal with a fractal di- 
mension df — 2, while for (7 = 1 the network is small- world 
and non-fractal with an infinite fractal dimension. Thus, 
when parameter q increases from to 1, the network fam- 
ily exhibits a crossover from "large" to small worlds, and 
simultaneously undergoes a transition from fractal to non- 
fractal scaling, i.e., its fractal dimension grows from 2 to 

CXD. 

The peculiar structure of the two network families makes 
it possible to explore the role of fractal dimension, instead 
of the power-law degree distribution, in the behavior of 
random walks on scale-free networks, since some particu- 
lar cases of them share the same degree distribution. This 
is the main goal of this paper. 



3.1 Average trapping time for (x, i/)-f lowers 

We study how the trapping time JT^"^ changes with n. 
According to the construction algorithm, for an edge in- 
cident an old node i and one of its old neighbors j at 
iteration n, it will generate x + y ~ 2 new nodes at it- 
eration n + I, among which x — 1 nodes, denoted by vi, 
V2, - ■ -, Vx-i, together with nodes i and j form a path 
{i, vi,V2, ■ - ■ , Vx-i,j) of X links, and the other y — 1 nodes, 
denoted by ui, U2, - ■ - , Uy^i, together with i and j con- 
stitute a path {i,ui,U2, ■ - ■ , Uy-i,j) of y links. We call vi 
and ui first-order (direct) neighbors of i, V2 and U2 second- 
order neighbor of i, V3 and U3 third-order neighbor of i, 
and so on. 

Let's examine a random walk on Fn+i{x,y)- By con- 
struction, upon growth of the networks from generation 
n to next generation n + 1, the degree of node i doubles, 
that is, it grows from ki{n) to ki{n + 1) — 2ki{n)- Among 
these 2 ki{n) neighbors, one half are the first-order neigh- 
bors of i with each belonging to a path of x links, and the 
other half are the first-order neighbors of i , that belong to 
the paths of y links. Let Z be the FPT going from node i 
to any of its ki{n) old neighbors, i.e., those nodes directly 
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connected to node i at iteration n; let X; (1 < ^ < 2: — 1) 
be the FPT originating at any of i's ki[n) Zth order new 
neighbors belonging to a path of x links to one of i's ki [n) 
old neighbors; and let Ym (1 < to < y — 1) be the FPT 
for going from any of the ki{n) mih order new neighbors 
of i, belonging to a path of y links, to one of i's ki{n) 
old neighbors. These FPTs obey the following backward 
equations: 



Z 

X2 



Yi 
Y2 



Yy^l 



-Z) + i(l+X2), 
Xi) + i(l + X3), 



= W 



1(1 



-id 

hi(i 



^ Xx-2 
1(1- 



Y2).. 



(5) 



Yv-2) . 



Equation ([S]) has a solution Z = xy. Thus, upon the 
growth of the (x, y)— flowers from generation n to genera- 
tion n + 1, the FPT from an arbitrary node i to another 
node j (both i and j are already existing in Fn{x,y)) in- 
creases by a factor oi xy. Then, we have = xyT^"'\ 
which is a basic feature of random walks taking place on 
the (x, 2/)-flowers and will be very important for the follow- 
ing derivation of the exact formula for the average trap- 
ping time. 

In order to determine the average trapping time (T)„, 

one may alternatively evaluate J^ien T^^\ To this end, 
we introduce some variables that will be convenient for 
description. Let /?„ denote the set of all nodes in Fn{x, y) 
that enter the networks at iteration n. Then, we can fur- 
ther define the following intermediate quantities for 1 < 
g <n: 



T, 



(n) 



5, tot 



and 



T, 



(n) 



5, tot 



(n) 



(6) 



(7) 



Hence, the problem of determining (T)„ is reduced to find- 
ing T^"of It is evident that f2„ = i7„ U f2n-i, which leads 
to 



rp{n) _ rp\ 

n.tot n — l,tot 



T. 



(n) 



-'n-l.tot 



(") 



(8) 



new nodes obey the following relations: 



J-vi 



(n+l) 



rTn(ri+l) 

rp{n+l) 
J-ui 



(n+l) 



T, 



rp(n+l) 




1 1 _L rp{n+l}^ 

2 1 r -r J 



M 1 + Ti. 



(n+1)^ 



(9) 



i(i + ri:+;))) + i(i + r("+^) 



From Eq. we have 

+ Tj,l'+p - 2{x - 1) + T^"+^^ + Tj"+^) , 
^ j,(„+i) ^ - 3) + Tif -I- Ti;:'+i^ ' 
+ = 2{x - 5) + + Ti^"_+/) , 

and so on. From these relations, we obtain 



1=1 



(„+l) ^ X{X^ - 1) ^ X - 1 



6 



(n+l) 



(10) 



T. 



(n+l) 



In a similar way, we can derive that 
^ 6 2 V * 

m— 1 

(11) 

Summing Eqs. (jlOp and ([Tl]) over all the M„ edges 
pre-existing in Fn{x, y) leads to 

^(n+i) _ xjx"^ - i) + y{y^ - 1) , , 

-'n+l, tot — g "^^n 

^ X + y -2 ^(^+1) 



+ E (^^(") 2 

x(a;2 - 1) + y(y2 - 1) 



6 



(x + y)' 



{x + y-2)Tl^Z^ +2{x + y-2)Tri: 



(n+l) 



+ 2"-i(x + y-2)r[';+^) + 



2"-Ma: + 2/-2)To\' 



(n+l) 



tot 



(12) 



where we have used the relation T'^*-"-' xyT^^^ ^\ Equa- 



tion 



indicates that before calculating T"^" ot , 



we should 



first evaluating the quantity r„ that can be derived as 
follows. 

As mentioned above, for each edge attaching i and j in 
Fn{x, y), it will generate x + y — 2 new nodes at iteration 
n+l, say wi, V2, ■ ■ ., fx-i, ui, U2, • ■ Uy^i. For a random 
walk in Fn+i{x, y), the trapping times of these x + y — 2 



Analogously, we have 

n(n+2) _ x{x^ - i) + y{y^ - 1) 



T. 



n+2,tot 



6 



(x + y) 



n+l 



{x + y^2) Ti'ltit + 2{x + y-2) f^^.t? 



(n+2) 



+ ... + 2"(a: + y-2)Ti\,,, 
2-{x + y-2)fl,;tt^. 



(13) 
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Equation ([T5| minus Eq. (fT^ times 2xy and consider- 
ing the relation — xyTj;^^^\ we obtain the recur- 
sion relation: 



Thus, by definition given by Eq. the average trapping 
time {T)n is 



^ ^-{x + y- 2xy){x + y)" ■ 

6 

(14) 

Using the initial condition T^^^^^ = ^^'^^'^ -i^+v) _ 
Eq. (ITil) is solved inductively to obtain 



(«+i) 



xy-l 



^ ^i"ot ^ Gi / a: + y - 1 _ x + y 
7V„-lVj; + y-2 " x + y -2) 

In[xy(a;4-y)] 

+ 7V„-lU + y-2 " x + y-2) ^'^^ 

In the limit of large network size, (T)„ is approximatively 
equal to 



+ C2 [xy{x+y)Y 



xy-l 



(15) 



(20) 



in which Ci and C2 are two constants with Ci — 
y - 2xy) and = f[%, = i^+vf-i^+v) _ ^y. 
Inserting Eq. into Eq. ([5]) yields 



ri",L = xyTt\'L + ( 3^ + C2 ) [xy(a; + y)] 



xy-l 



n-l 



Ci 



xy-1 



{x + yY 



(16) 

Considering T|Vot = ^''^yi p^'+yi ^ Eq. (HH) is resolved by 
induction to obtain 



^/Ecpiation ([201) shows that, the dominating scaling of 
average trapping time grows as a power-law function of 
network size Nn with the exponent 9{x,y) = Inxy/ ln(x + 
y) being compatible with the bounds derived in |31| . Note 
that a scaling similar to 6{x, y) was previously obtained for 
general fractals without the scale- free phenomenon [T^ . It 
is also worth emphasizing that for some limiting x and y, 
d{x, y) recovers some previous results. For x = 1 and y = 
2, 0{x,y) reduces to In2/ln3 previously obtained in |27) : 
for a; = 1 and y — 3, 0{x,y) is reduced to In3/ln4 |48) : 
while for x — 2 and y — 2, 9{x,y) corresponds to 1 [48) . 
The consistency confirms that Eq. (PO)) is valid. 



Result analysis 



rp{n) _ 
n,tot 



Do 



-(x + y)"+ 



xy ~ X — y 
D3 D2{x + y) Di{x + y) 



xy xy{xy - x - y) xy{x + y - 1) 

+ , f ' 7^[xy{x + y)T ■ 

xy{x + y-l) 



where i^i = ^ + C2, D2 = and = T[ 

{x+y)^-{x+y) 



(17) 

(1) _ 

tot ~ 



From Eq. we have {x + y)" = fqrf^^n 



x+y-2 



and n = log™,,, ( i Nn ■ In addition, we as- 

^x+y \x-\-y — 2 x+y — 2 J ' 

sume that Gi = G2 - ^ ^^^^ £i(£±£L 



xy—x—y ' 



xy xy{xy — x — y) a;y(a;+y — 1) ' 

and G3 = -zr(-^7~zr\- Then, Eq. ffTTll can be recast as a 

xy[x-\-y 1) 

function of network order 



x + y — 2 x + y — 2 



ln[a;-y(3;-f y)] 

.x + y — 2 x + y — 2/ 



As mentioned above, the average trapping time (T)„ is 
a highly desirable quantity characterizing the trapping 
problem in a network. For trapping issue on the (x, y)— flowers, 
(T)„ is dominated by the exponent 0{x, y) = It is 

not difficult to verify that for different values of x + y, the 
exponent 9{x, y) increases with the sum of x and y. Hence, 
x + y determines the trapping efficiency for random walks 
on the (x, y)— flowers: the larger the sum x + y, the less 
the efficiency. This can be understood from the following 
arguments. When x + y become larger, the (x, y)— flowers 
are more homogeneous, since the exponent 7 of the power- 
law degree distribution grows with increasing x + y. Ac- 
cording to the previous result that the trapping efficiency 
in scale-free networks decreases when the exponent 7 in- 
creases [26j . we can know why the sum x + y influences 
strongly the leading behavior of the average trapping time. 

For the case when x + y is flxed, Eq. (|20p shows that 
(r)„ is determined by the difference between x and y 
(since we have assumed that x < y) : the smaller the differ- 
ence y—x, the larger the value (T)„. As shown in section[5J 
for the (x, y)-flowers with fixed x + y, they have the same 
degree sequence and thus the same degree distribution. 
Thus, we can conclude that power-law degree distribution 
alone is insufficient to determine the trapping efficiency in 
scale-free networks. 

As pointed out above, the power-law degree distri- 
bution alone cannot determine the behavior of trapping 
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problem in scale-free networks. Then, it is natural and 
important to ask: which structural property plays a cru- 
cial role in the trapping efficiency in the (x, 2/)-flowers with 
the same value of x + y7 Since for those {x, y)-flowers with 
an identical sum of x + they have the same average 
node degree, the same degree distribution, and the same 
(zero) clustering coefficient; moreover, for the family of 
(a:, y)-flowers at the same generation, their initial nodes 
have an identical degree, we argue that the difference for 
the average trapping time shown in Eq. (j20[) is only due 
to the fractality. For those {x, ?/)— flowers having the same 
parameter of x + y, when x increases from 1 to its maxi- 
mum, their fractal dimension df — '"j^^^'' decreases from 
00 to its minimum, while the exponent 9(x,y) increases 
from its minimum to maximum. Concretely, the trapping 
efficiency in the {x, y)— flowers with the same x + y is dom- 
inated by the fractal dimension: the larger the fractal di- 
mension, the more efficient the trapping process. Partic- 
ularly, when x = 1, the network has an infinite fractal 
dimension, which corresponds to the most efficient struc- 
ture for trapping, among all (x, ?/)— flowers with the same 
sum of X and y. 

It is worthy of mentioning that the impact of fractal 
dimension on the trapping efficiency can also be seen by 
comparing the average trapping time in the (1, 2)— flower |27j 
and the two-dimensional Sierpinski gasket . The trap- 
ping process is more efficient in the (1, 2)— flower than in 
the two-dimensional Sierpinski gasket, in spite of the fact 
that both networks have the same numbers of nodes and 
edges at any iteration. The difference for average trapping 
time in the two networks is at least partially attributed to 
the fractality. The (1, 2)~flower is non-fractal with an infi- 
nite fractal dimension, while the corresponding Sierpinski 
gasket has a fractal dimension of 1 + 



that was generated by the first iterative method and is 
attached to this new neighbor of i, to an old neighbor of i; 
and let D be the FPT originating at any of i's new neigh- 
bors created by the second iterative method to one of i's 
old neighbors. Then the following backward equations can 
be established: 



C = 1 + 1(1 + 

1 + 1(1 + A), 



(21) 



which leads to A 



relation t/"+^^ = J2 

q+3 



^q^. Thus, we can obtain an important 



Having obtaining the evolution rule of trapping time 
when the hybrid networks grow, we now compute the the 
average trapping time (T')„. For simplicity. In the follow- 
ing text, we will use the same notations as those for the 
{x, y)— flowers deflned above. Analogous to (x, y)— flowers, 
we have the following equation: 



rp{n) _ rp\ 

n.tot n — l,tot 



T. 



in) 



12 



2 n— l.tot 



(71) 

Hence, the problem of determining T„ is reduced to 

evaluating T^'tot that can be obtained as follow. 

According to the iterative rule, at a given generation, 
for each edge connecting two nodes i and j in iJ„, it will 
generate two new nodes at the subsequent iteration n + 
1. Both nodes are either generated by the first iterative 
method with probability q and labeled by wi and W2, or 
generated by the second iterative method with probability 
1 — q and denoted by w[ and w'2- For wi and W2, their 
trapping times obey the relations: 



(") 



(22) 



3.2 Average trapping time for hybrid networks 



(«+!)_! 



(«+!)_ 1 



("+i)i 



W2 J 
("+l)l 



^[1 
^[1 



(n+1)] 
("+l)l 



(23) 



We proceed to show that the above relation between av- 
erage trapping time and fractal dimension for trapping 
problem in the deterministic scale-free (a;, ?/)-flowers are 
also true for the random hybrid networks 

Similar to the {x, ?/)-flowers, we first establish the re- 
lation dominating the evolution for T^^"^ with iteration 
n. Let's consider a node i in iJ„. By construction, we 
know that, upon the growth of the networks to itera- 
tion n + 1, the degree of node i increases from ki{n) to 
ki{n -I- 1) = 2ki{n). Among these ki{n + 1) neighbors, 
some are created at iteration n + 1, each of which is either 
generated by the first iterative method with probability 
or created by the second iterative method with probabil- 
ity 1 — q. We now consider an unbiased random walk in 
Hn+i- let A denote the FPT for going from node i to any 
of its kiin) old neighbors, i.e., those nodes directly con- 
nected to i at iteration n; let B stand for the FPT from 
a new neighbor of i, which is generated through the first 
iterative method, to one of its ki{n) old neighbors; let C 
represent the FPT from a new node created at iteration 
n + 1, emerging simultaneously with a new neighbor of i 



From Eq. ([23)) . we obtain 

^(n+l) , rp(n+l) _ 4 I rp(n+l) rp(n+l) 



(24) 



Similarly, for w'l and w'2, we have 

T^?+i) = t/^?+i) = 1 + 1 + 7;("+i)) , (25) 
yielding 

rp{n+l) _^ y(n+l) _ 2 + 4. rp{n+l) ^^q^ 

Summing Eq. and Eq. over all the M„ edges 
in Hn and considering the iterating probability for each 
preexisting edge at the generation n, we get 



filt% = {2q + 2)A/„ + (h{n) X T, 



-i(n+l) 



{2q + 2)4- + 2Tt+J^+2%t+^l, 



.2«T,(«+i)+2"To("+^). (27) 
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In addition, from Eq. (P7)) . we can write out T^+t^tot 

AU + l 



(n+2) r,2T=.(n+2) 
ri+l,tot + ^ -'ri.tot 



_9ri+lrf.(n+2) r,n+lr^=,(n+2) 
^ -'I, tot ^ -'O.tot 



(28) 



Equation (|^ minus Eq. ([27]) times and apply- 



ing the relation of 7^^("+^' 
following recursive relation: 



q+3 

-i^T.'"+^\ we obtain the 



^n+2,tot ~ „ I q^n+l,tot - ^ ^ 



q + 3 



q + 3 



(29) 



With the initial value of T^l^^^ = 6 + q, Eq. ^ is solved 
inductively 



(") 



48 + g2 _ ^ / 48 



(g + l)(6-2g)^„_i 



9-g V3 + g/ 9-g 
Plugging Eq. ^ into Eq. ([22]) yields 

.(n) _ ^2_ („_i) ^ 48 + - g / 48 

-'n-l,tot 



(30) 



q + 3 



9-q 



3 + q 



(g + l)(6-2g) 
9-q 



(31) 



Considering T^^l^ = 10, Eq. §^ can be resolved by in- 
duction. When 5 = 0, 



.(«) 



18' 



x4" + 3n+10); 



(32) 



when < g < 1, 

,(„) _ (g + 3)(48 + g^-g) 



T 



48 



36(9 -g) V9 + 3 

(3-g)(g + 3)(g+l)^„ 



+ 



2q(9 - q) 
(g + 3)(g^-g + 3) 
ISq 



12 

3 + g 



(33) 



Using Eq. Q, we can find the analytical solution for 
average trapping time in _ff„, from which we can obtain 
the exact dependence relation of {T)n on the network size 
Nn. For q^O, 



{T)r 



3iV„ - 4 



36(7V„ - 1) 
while for < g < 1 



127V„ 



3 1n(|jV„-2) 
2 In 2 



6 



; (34) 



gT3 



m 1 (g + 3)(48 + g^-g) ^3^^ _ ^ - 

^^""iV„-l 36(9 -g) U 



1 (3-g)(g-f 3)(g+l) /3 



7V„ - 1 2g(9 ~q) V 2 

1 {q + 3){q^-q + 3) f3 



Nn-2 



Nn-1 



18q 



(35) 



1.00 



0.95 



0.90 



3 



0.85 




0.0 0.2 0.4 0.6 0.8 1.0 

q 

Fig. 5. (Color online) Relation between 6{q) and q. 



From Eqs. and we can easily see that in the 
full range of < q < 1, the average trapping time (T)„ 
in the hybrid networks behaves as a power-law function of 
network order Nn, with the exponent 9{q) = ln^q^/ln4 
decreasing with parameter q. Concretely, when q grows 
from to 1, 9{q) drops from 1 to In3/ln4, see Fig. [5] 
In addition, when q enhances from to 1, the fractal di- 
mension d{ of the hybrid networks increases from 2 to oo. 
Thus, in the random hybrid networks, the diffusion effi- 
ciency also enhances with the increasing fractal dimension 
as that found for the (a;, y)— flowers. 



4 conclusions 



We have studied the trapping problem in two families 
of networks, i.e., the (a;, y)— flowers and random hybrid 
networks. The former is deterministic, while the latter is 
stochastic. Both display rich and specific structural prop- 
erties and can recover some prominent features of various 
real networked systems. We focused on a particular case 
of trapping issue with the immobile trap located at either 
of the initial nodes of the networks. We derived analyt- 
ically the average trapping time and demonstrated that 
they grow as a power-law function of the network size 
with the exponent varying with network parameters. We 
showed that the power-law degree distribution (even de- 
gree sequence) itself is not sufficient to determine the be- 
havior of average trapping time of random walks on fractal 
scale-free networks. Instead, we presented that fractal di- 
mension plays a predominant role in the scaling of average 
trapping in both families of fractal scale- free networks un- 
der consideration. Finally, it should be stressed that here 
we only studied particular networks with the trap fixed 
on a specific node having the highest degree, whether the 
conclusion also holds for a general fractal scale-free net- 
work or for a non-hub node needs further research in the 
future. 
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